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Abstract 

The relation between the ratio T ee (ip )/T ee (J/ip) and T 27 (?7 c )/r2 7 (?7c), expressed in 
terms of the configuration mixing amplitudes induced by the contact spin-spin in- 
teraction of quarks in the ground and radial excitation states, is shown to give, after 
the inclusion of the newly derived relativistic corrections, the radiative t] c - and 
?7 c -width ratio in fair accord with recent experiments. The dynamical model is pro- 
posed to derive the ratio of relative probability of the ground (r/ c (2980))- and first 
radial excitation (?7 c (3640))-state formation in 77-collisions followed by their decay 
into the KKir and pp channels. 



1. Two-photon decays of heavy quarkonia provide valuable information on heavy 
quark dynamics and have been under consideration in recent experimental PU El El IH El El 
El El E] and theoretical (121 UB UH UB E3 HZ1 HE] studies. It seems reasonable to assume, 
at least as a first approximation, that there is no mixing between light and heavy quark 
sectors, and that one can consider the needed transition amplitudes separately for mesons 
constructed of the u-,d-, s-, and heavy c- and b-quarks. 

In this work, we concentrate on the charmed quark sector and our main concern in 
this problem will be the question of the degree of model (in) depen- 
dence of the S-wave QQ annihilation rates, or their ratios, with respect to the role of 
short-range spin-dependent forces. This is an important theoretical question because 
in many approximate relativistic approaches to description of the annihilation of bound 
antiquark-quark S'-wave states that are subjected to strong short-range interactions one 
needs to introduce the cut-off procedures to get rid of singular behaviour of matrix ele- 
ments considered, the cut-off parameters or the "smearing" procedures being introduced 
basically on the phenomenological grounds. 

2. We start with just postulating for a mass operator of the heavy quark (Q = c, b) S'- 
wave systems the simplest spin-dependent matrix 

(n\M\n) = MJ nn , + C ■ 1^(0)^(0) ■ (ff Q ■ S Q ) (1) 

and write down simple perturbation theory expressions to have an estimation for the 
mutual change of the ground and first radial excited states due to switching on the contact 
spin-spin potential 

MO) = MO) + E \m)(m\V\n)/(EV> - E^), (2) 
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where \i) = ipi, and are the wave functions and energies of the unperturbed hamil- 
tonian and V is the perturbation potential that we identify with the contact spin-spin 
interaction of quarks. 

Applying Eq.(J2J) successively to ipns(J), where n — 1,2- are the quantum numbers of 
the S-wave radial excitations of the cc- quarkonia with the spin J = 0, 1, using the spe- 
cific (i.e., proportional to the ^-function form of the perturbation potential) and keeping 
everywhere the terms of the first order in 0(V) we get the relations 



R 2S/ls\hfs - , T .J f n \2\ h f» 



fe(0) 2 [(1 + V^ml - mg)- 1 + V£(m° 2 - m°)-' + . . .)] 2 
<M0) 2 ' [(1 + V£{ml - m")- 1 + V&(m\ - m°)-l + . . .)] 5 
fe(0) 2 



[1 + 2(s J 2 - s(], (3) 



where s^, J = 0, 1 is the sum of terms of the order 0(V) entering into the denominator 
and numerator of (3). Using the evident relation V J=0 = — 3V J=1 we obtain 



R 2S/ls\hfs + 3 • R 2 s/is\hfs — 4 ' ^ 15 (0) 2 ' ^ 

We make an estimation for the ratio T(rj' c — > 77) /T(r] c — > 77) on the basis of experimental 
data of leptonic charmonium decays, approximate validity of the lowest order perturbation 
theory for the color-hyperfine splitting interaction as well as on our anew obtained form 
of the relativistic corrections to the considered decays. 

We remind first the known results PUJ for the lowest order QCD corrections 



ml T(n c — > 77) 4, a s , ,,,, 

' - (1 + 1.96— ) 1 VcK n 2 , (5) 



m j/v r (W^ e+e ) 3 71 \^j/i> 

where a s should be evaluated at the charm scale [TEj . 
In certain chosen ratios such as 

r & r ^ = ( w ? rJ 2s?is (6) 

V^Tt'e m v ,m J/lP Rig/is 

the ratios of ip n s(OYs are cancelled and also the QCD radiative corrections are assumed 
to be mutually compensated not only in the next-to-leading order (NLO), but also in the 
higher orders of the perturbation theory, which, as seen from refs. |15tlT^ I20j. are not neg- 
ligible. We assume further on that each type of studied corrections can be approximately 
represented in factorized form 

*ns(0) = hMO) | 2 (1 + 5\rad) + 5 J nS (rel) + 5 J nS (hfs)) ~ 

~ |^ 5 (0)| 2 (1 + 6 J (rad))(l + 6 J nS (rel))(l + 6 J nS (hfs)). (7) 

That means that the 1 + 5 J (rad) factor coincides in the lowest order radiative correction 
with 1 + (a J /ir)a s (m c ), where (a J=1 = —5.34) and (a J=0 = —3.38), and it is assumed to 
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be cancelled, together with higher order terms, in both the leptonic {R^shs) anc ^ 2-photon 
(Risns) "single ratios". 

The notation S^ s (hfs) refers to the correction factor due to the spin-spin potential- 
induced factors in (jH) and this type of corrections is not cancelled in the double ratio (jHJ) 
unlike the mentioned linear combination of the "single ratios" 

The relativistic correction factor 1 + S^ s (rel) is defined in the following manner. The 
static approximation for both T ee and T 77 , resulting in their proportionality to the respec- 
tive wave function value "at origin", ip n s{0), follows from neglecting of the dependence of 
the bound quark annihilation amplitudes on their internal motion momenta 

| J dpA(cc(p) - ee( 7 7)) <Pns(p)? ~ \^ n s(0)\ 2 \Athr(cc - ee( 77 )| 2 . (8) 

We have found that the appearance of the ip n s{ty together with the relativistic correc- 
tion factor 1 + 5^ s (rel) can be justified from more general expression that includes the 
non-static annihilation amplitudes depending on the internal quark momenta taken into 
account. 

2^=M) * g/^„ sW ^ log (M±£)// l^fr), (9) 

v /I ^^-/|W p)(1 -^ )/ /|W rt ' (10) 

where <p n s{p) is the nS-state wave function in the momentum representation, and e n s is 
the quark energy in the nS'-state. The momentum- dependent factors of the (cc) n s —* 
e + e _ ( 77 ) amplitudes are given, e.g., in [T7], but instead of taking e(p) = (m 2 +p 2 ) 1 / 2 we 
prefer to define the continuation of the cc -annihilation amplitudes to the bound state 
kinematics following the so-called "on-energy-shell / off-mass-shell" prescription when 
e(p) — > e n s = m n s/2 now remains independent of the internal motion momentum \p\ 
while being averaged with the wave functions <p n s{p) i n (fTUj) . This picture of the bound 
state dynamics underlies the derivation of the relativistic Schrodinger-type wave equa- 
tions such as the quasipotential equation suggested by Todorov or different variants 
thereof, e.g., [21]. In this prescription, the relativistic correction factor for the electron- 
positron annihilation of vector charmonia is especially simple because there are no addi- 
tional momentum-dependent factors in the integrand, while 5%s°{rel) is derived directly 
from the relation 

oh , jrJ=0/ ;mV2 m c t d3 P a i \ e «s, [tnS + P,, f d 3 p 

JJl^ POO 

= 2— / dr^ nS (r)sin(e nS r) /ip nS {0) (11) 

CnS JO 
TY> 2 1 

z> 2^ " fe(0) - —^ 5 (0) + • • .]Mrf(0), (12) 

e nS € nS 

where masses m n s correspond to masses calculated without spin-dependent corrections: 
mis - (3mj/^ + m Vc )/4,,m 2 s ^ (3ny + "\/)/4. 

Keeping only the first term in asymptotic series for the Fourier integral (|12p. we obtain, 
in the accord with the Riemann-Lebesgue lemma, the simple approximate relation 

[i + ^°M)] 1/2 -^ (13) 

e nS 
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demonstrating more strong dependence on the relativistic corrections of the two-photon 
decay amplitude as compared to the leptonic one 



Besides the relativistic corrections to be taken explicitly, we have observed a signifi- 
cant contribution to our relation for ratios (jOJ due to inclusion of the short-ranged spin- 
dependent interaction which modifies the vector and pseudoscalar wave functions "at 
zero" quite asymmetrically. However, for the linear combination (@| of two ratios the 
"hyperfine" corrections are compensated up to terms of the order 0(V 2 ). 

Collecting now all found corrections we get the resulting relation between the widths 
of the lowest lying states of charmonia 

m v ,m 2 2S g r£ m^m 2S (m ls + m c ) = ^g(O) 2 
Hence, it follows that 

-g = 0.21 ± 0.06, (16) 

if we take (V>2s(0)Ms(0)) 2 = 0.653 and m c = 1.48 GeV according to j2HI221, , masses 
and leptonic widths from PQ. A rather large uncertainty of the ratio obtained is largely 
due to experimental errors of the measured leptonic widths. It should be noted that unlike 
the very ^> n s(0)'s their calculated ratios are much less model-dependent . In particular, 
the ratio entering into (fTo^) calculated with the Cornell (i.e., the " linear + Coulomb") type 
potential is close to that calculated with the "running" a s (r), which makes the one-gluon- 
exchange potential softer at small distances, although the very Vis(0)' s differs about 
two times from each other (221 • We believe that mild smearing or regularization of the 
short-range quark interactions providing formal finiteness of ip n s{0)'s following from the 
relativistic equations will also leave their ratios relatively intact. 

As has been mentioned, the ratios of ip n s{0) are dropped in © and one can obtain a 
kind of the lower bound for the double width's ratio leaving only one, beyond and next 
to unity, term in every infinite sum of the perturbation corrections to ^5(0) due to the 
contact spin-spin potential of charmed quarks treated as a first order term of perturbation 
theory. Including then already fixed relativistic corrections and assuming, as earlier, the 
cancellation of the (static!) radiative corrections, we obtain 

T^lte m rh m ls m^ (m is + m c ) 

yj=0 , yJ=Q V J=1 -\- V J=1 V J=1 -I- V J=1 

(1 + V_n +V22 )2/(1 + V_n +V22 )2 ^ 1 _ 8 ^ii +V22 
m 2S - mis m 2 s - mis m 2S - m is 



from where a new constraint follows 



r 



XL 

77 



>0.1. (18) 
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and where the numerical values PQ for mfg 1 , m^g 1 , defined earlier, and V 7 ^ -1 = {rrtj/^ — 
V 22 =1 = (m^i —m^ )/4:, were used. The comparison of (fTBj) and (JTSJl tells about 

a significant role of sums < nS\V J \nS > over n > 1, 2 in the definition of the individual 

ratios of T^/r£ and especially of Ffa/T*. 

To make contact with the available experimental data for radiative widths of charmonia 
it is necessary to estimate also their relative branchings referring to the studied hadronic 
decay channels. 

The main assumptions underlying our estimations of the branching ratios and 
B^ c entering into the experimentally measured processes of the two-photon fusion pro- 
ducing rj c and rj c and subsequent hadronic decays r\ c (rj c ) — > h , where h = KKtt [3] or 
h = pp are the following. We assume a simple kinematic structure of the respective 
decay amplitudes and an approximate dynamical assumption for the ratio of relevant 
couplings or, rather, complex form-factors in the considered vertices 

A{ Vc { V ' c ) -> KKtt) = g{ Vc { V ' c ) -> KK^-^rifiK^, (19) 

MVciv'c) PP) = g(Vc{v'c) PP)Fc(rn 2 c5 )u(P p )^ 5 v(P p ), (20) 

A(Vc(Vc) -> G m G rn ) = g(r) c (r) c ) -> G m G m )e tlvp(T e^qle l 2q2 , (21) 

Ifffa^ ^)| 2 _ ^)| 2 ^2) 

|c/(r/ c -> G m G m )\ 2 \g{v' c ^ G m G m )\ 2 ' 

where only isospin structure of the KKtt decay channel is indicated in the first line and 
the generalization to the SU (3)-symmetry can easily be written down. We include the 
form-factor F C (Q 2 ) in ({20} to mention about its possible variation depending on the time- 
like momentum transferred in the interval m 2 < Q 2 < m 2 , . We note that unlike the 

vector charmonia J /if)— and -^-decays, the pseudoscalar decay branchings B(r) c (r] c ) — > 
(h, 7)(cc)) containing the charmed quarks in the final state, e.g., the decay r\ c — > lTTr\ c , 
are much less significant due to smallness of the coupling constant a s (m c ) ~ 0.3, as 
compared to the branching ratios B(t] c (t] c ) — > h), where h denotes hadron states composed 
of the light (it, d, s) quarks. Following a usual practice (or, alternatively, the quark- 
hadron duality hypothesis for inclusive processes), we identify the total width of these 
processes with the width of the bound cc-quark annihilation to pair of "free" gluons 
F(Vc(Vc) ~^ G m G m ). Further, we attribute to gluons finite "effective" (or dynamical) 
mass mo of the order 0.7 GeV, which was advocated in ref. |2Sj arid in some earlier works 
cited therein, on the basis of the detailed study of the experimental photon spectrum in 
the inclusive reaction J /if} — > ^G m G m — > jX. The evaluation of the important transition 
probabilities T) c (Vc) ~* KKtt, studied in several experiments, see, e.g., [3] and further 
references therein, has been performed with the help of the integral relation for invariant 
3-body phase space 

p(m ^m 1 + m 2 + m 3 ) = 1 2 f —F(s, s u s 2 , s 3 , s 4 ), 

F(s, si, s 2 , s 3 , S4) = [{s - si)(s - s 2 )(s 3 - s)(s 4 - s)] - 5 (23) 

where s 2; i = (mi ± m 2 ) 2 , = (m ± m 3 ) 2 and we choose m = m Vc (r l ')^ ' m i = 

m 2 = w,Ki ^3 = "V- After that, the integral which is the essential part of the ratio 

$kk I^kk can easn y be calculated numerically. 
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Including standard relativistic normalization factors of the initial states, summing 
and averaging over spin degrees of freedom of the listed amplitudes squared and using the 
assumed relation (J22j) . we obtain 

&l° 0.83, for h = KKtt, 

Bf ~ 0.65 ■ (F c (my /F c (m 2 VB )) 2 , for h = pp. (24) 

Hence, our results are presented in Table 1, where the upper bound for the pp -decay 
channel refers to the ratio of the form-factors put equal to unity and where we have 
included also some recent theoretical predictions and experimental results for the widths 
of the Vc(Vc)~^°~11 decays and their ratios. 

Table 1: Recent theoretical and experimental results of the rj c and r] c two-photon decay 
width (the notations: h = KKir for ref.jH] and h = pp for ref.JTJ E] , the entry with su- 
perscript (x) is our estimate using other works) 





KeV 


T% KeV 


1 27/ 1 27 


1 / J- 27°^ 


PDG [T] 


7.5 ±0.8 








CLEO 12 E] 


7.6(0.8)(2.3) 






0.18(0.05)(0.02) 


E760/E835 E] 


6.71^(2.3) 






<0.16 


L3 


6.9(1. 7)(2.1) 


< 2 






DELPHI 






< 0.34 




Gupta [12! 


10.94 








Miinz[T3] 


3.50±0.40 


1.38±0.30 


0.39±0.10^ 




Chao dU 


6-7 


2 


0.28 - 0.33^ 




Fabiano < 


S.18(0.57)(0.04) 








Ebert |17j 


5.5 


1.8 


0.33^ 




Kim|IBl 


7.14±0.95 


4.44±0.48 


0.62±0.10^ 




This work 




1.6±0.5 {a;) 


0.21±0.06 


0.18 ,h = KKix 



< 0.15, h = pp 



The behaviour of the unitary- singlet, pseudoscalar form-factors in the hadronic tran- 
sition vertices is of considerable interest for the understanding of mechanisms of se- 
quential processes rj c (7] c ) — ► G m G m — > light hadrons. The closeness of our estimated 
R(r) c /r] c ) decay-ratio in the KKir channel, with no additional form-factors included, to 
the CLEO data looks intriguing and needs more investigation to be understood. At any 
rate, in the pp channel one should expect a stronger dependence of the result on the 
(m 2 c /m 2 , ) n ratio, where the effective power n is expected to be two units larger in the 
decay amplitude for two-baryon final state as compared to the transitions into the two- 
meson states, according the quark-counting rules. The two-gluon state mediating the cc - 
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and q n q n - states, where q = u,d, s, would provide a new testing ground for checking the 
generalized parton approach or the diquark model which were successful in the description 
of the two-photon annihilation processes, like 77 — > pp, etc. Therefore, a further study 
of the reactions pp — > T] c {Vc) ~^ 77 ( or KKt:) with better statistics and accuracy is of 
interest, e.g., via the ongoing experiments at FNAL or at a planned antiproton storage 
ring at GSI. 

3. Finally, we note that relation (|15|) can be applied to any pairs of the "hyperfine- 
split" radial-excited states. In particular, using it for pairs of ratios R^s/is an d R3S/2S with 
the needed input values = 4039 MeV, r£ = .89±.08 keV |2Z1I2H|, |^ 3 s(0)/^i5(2S)(0)| 2 
= .56(.86) 122] we obtain the estimation ~ 4003 MeV and T^/T^ ~ .22 of the 
"naked" (i.e., without possible hadronic corrections due to virtual, open-charm interme- 
diate (DD* + DD*) - states) parameters of the still to be observed 77^ -resonance. 

Our main results, relation (j!5|) and the numerical entries in Table 1, demonstrate 
a considerable suppressing effect of the relativistic and "hyperfine" spin-dependent cor- 
rections on the recently observed two-photon decay of the r/ c (3640) -resonance. If it 

is true, this effect should also display itself in the total width Y^ t represented by the 
decay into two gluons, either massless or effectively massive. Using the average value 
Yf ot = 32.3 ± 2.2 MeV of the CLEO and BaBar g] results, we have got the estimate 
of the total width of the r/ c (3640)- resonance 

(32.3 ± 2.2) ■ (.21 ± .06) ~ 6.8 ± 2.0 MeV, for m G = 0, 
r£ = (32.3 ± 2.2) • (.21 ± .06) • ((1 - 4mh/m^)/(l - 4rr4/r<)) 3 / 2 
~ 7.8 ± 2.3 MeV, for m G ~ .7 GeV 

A more precise measurement of this important parameter, as compared to the published [Sj 
result = 6.3^4 2 4 MeV, would be very desirable. 
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